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Abstract. We show that all long nulls in PSR BOS 18- 13 are trains of rapidly alternating nulls and pulses (each shorter than 
one pulse period). Sometimes only the nulls coincide with our pulse window, resulting in one of the apparently long nulls seen 
occasionally. We show these are seen as often as expected if during such a train the probability for nulls is 1.2 times less than 
for pulses. 

During nulls, the subpulse drift-speed appears to increase. We assume that the carousel of sparks that possibly underlies the 
subpulses actually slows down, as it does in similar pulsars like PSR B0809+74, and conclude that the subpulse-drift in this 
pulsar must be aliased. The carousel must then rotate in 30 seconds or less, making it the fastest found to date. 
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1. Introduction 

Observations of drifting subpulses are a nice example of under- 
sampling: in some radio pulsars the positions of the pulse com- 
ponents change regularly (as in fig. ^, but as these positions 
are determined only once per pulse period, the exact underly- 
ing motion remains unknown. 

Finding this underlying motion, and its relation to other 
pulsar parameters, could be helpful in determining which, if 
any, of the proposed mechanisms to generate dr ifting subpulses 
JRuderman & Sutherland! 1 19751 IWrightll2003l) is correct. We 
will assume the subpulses are formed by discrete locations 
of emission ('sparks') that rotate around the magnetic pole 
of the pulsar, causing the pulse components to move through 
the pulse windows ('drifting subpulses'). Different directions 
and speeds of this rotation can produce identical subpulse po- 
sitions at the observer, so-called aliasing (see also fig. 2 in 
Ivan Leeuwen et alJ2003h . To determine the actual speed of the 
spark carousel, one can check for the periodicity caused as 
brighter or offset sparks re-appear after one carousel rotation. 
[Oeshpande & Rankin ( 1999) did so for PSR B0943+10, and 
found a carousel rotation ti me of around 41 secon ds, longer 
than expected from theory wuderman & Sutherland 1975). 
With the rotation time known, Deshpande & Rankin ( 1999) 
could track individual sparks. They noted that sparks retain 
their characteristic brightness and positions for about 100 sec- 
onds. If this number were the same in other pulsars, carousel 
rotation times longer than 100 seconds could not be found by 
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a periodicity search, as the characteristic traits of the subpulses 
will have changed before they return into view. For such long 
carousel rotation periods a different technique is needed. 

In some pulsars there is a mechanism that can be used to cir- 
cumvent the undersampling problem and hence determine the 
underlying carousel-speed: nulling. In these pulsars, a few per- 
cent of the pulses are missing as the pulsars suddenly turns off 
('nulls'). After these nulls, the subpulse drift is often affected. 
Van Leeuwen et al. (2002, 2003) investigated this drifting- 
nulling interaction in PSR B0809-I-74 and find it can only be 
explained if the underlying carousel rotates in over 200 sec- 
onds, again much longer than expected. They predict that in 
pulsars with shorter carousel rotation times the drift direction 
should reverse after nulls. 

In their 1983 paper, Lyne & Ashworth mention the nulling- 
drifting interaction of PSR B0818-13 and show the post-null 
drift behaviour (their fig. 14). For long nulls, there appears to 
be a reversed drift direction. Combining this with the prediction 
of van Leeuwen et al. (2003) we decided to further investigate 
the drifting-nulling interaction in PSR B0818-13. 



2. Observations and data reduction 

We have observed PSR B0818 -13 from 2000 t o 2003 
with PuMa, the Pulsar Machine JVoute et al.1 l2002h at the 
Westerbork Synthesis Radio Telescope (WSRT). The 1 1 obser- 
vations amount to 20 hours in total or, using the 1.238 s period, 
5.8x10"^ individual pulses. All observations were conducted at 
a time resolution of 0.4096 ms and around frequencies in be- 
tween 328 and 382 MHz, with a bandwidth of several times 
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Fig. 1. Stacked single pulses with two 
intermittent nulls in lighter grey. 



10 MHz, depending on interference. Each 10 MHz band was 
spUt in 128 channels and dedispersed offline. 

The data were searched for nulls and subpulses in the way 
described in van Leeuwen et al. (2002): for each observation 
we compose a histogram of the energy of individual pulses, like 
fig- 12 The discrete population around zero-energy pulses is la- 
beled 'nulls'. In each individual pulse, we fit Gaussian profiles 
to the subpulses. We determine which subpulse-fits are signifi- 
cant and store the subpulse positions, widths and heights along 
with the raw pulse data. As the subpulse-drift bands form regu- 
lar structures, we use the predictive capability of the thus-found 
strong subpulses to detect weaker subpulses over the noise. 



train to be the start and end of the group. In fig. 0] we give an 
overview of all event lengths. From this figure it is clear that 
the occuiTence of null-groups as a whole is not governed by 
chance: there are too many long events. Within an event the 
pulsar blinks fairly quickly, but the change back to the normal 
pulse state usually takes several seconds. 

The null group in the left panel of fig. [2 for example, we 
classify as an event of length 5, while the one in the right panel 
is an event of length 6. As the first and last nulls act as markers, 
there are only 6-2=4 positions in the rightmost of fig.[2group 
that can be either nulls or pulses. In this case, there are two 



3. Results: analysis of nulls and subpulse drift 

Statistics of short and long nulls 

PSR B0818-13 nulls about 1% of the time. As the histogram 
of null lengths fig.|3lshows, most of these nulls are very short, 
with lengths of only one or two periods. What the histogram 
doesn't show however, is that although roughly 2/3 of these 
short nulls occur by themselves, many are observed in groups 
of intermittent pulses and nulls, li ke the groups seen in fig . 
\l\ As this behaviour is uncommon, iLvne & AshworthI lll983t) 
chose to disregard all intermittent short nulls, and only inves- 
tigate longer nulls. But has the pulsar really turned off during 
such a long null? Or are the long nulls in PSR B0818-13 actu- 
ally intermittent nulls where the series of short nulls happen to 
line up with our line of sight? This is the question we investi- 
gate. 

If we assume that, during a train, the pulsar quickly and 
successively turns on and off", we can use the statistics of the 
null/pulse occurrence within a train to determine the likeli- 
ness ratio of null and pulse occuiTence. To classify the differ- 
ent pulse trains, we have defined the first and last nulls of a 
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Fig. 2. Energy histogram for PSR B0818-13. In Hght grey 
we show the energy distribution of the energies found in the 
on-pulse window. In black we show the energies found in an 
equally large off'-pulse window, caused by system noise. The 
set of low-energy pulses detached from both the normal and 
the null distribution are probably half-nulls. 
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Fig. 3. Histogram of the null lengths. 




Event length (pulses) 
Fig. 4. Histogram of the event lengths. 



nulls and two pulses in this middle part, in a NNPNPN pattern. 
In table [l] we give an overview of all the different null-pulse- 
trains present in our sample, in the same pattern notation. In 
the middle parts of all events, there are 238 possible locations 
to be divided between pulses and nulls. Of these, 107 are nulls, 
and 131 are pulses, so we conclude that, on average, pulses 
are 1.2 times more likely to occur within a group than nulls. 
Using this ratio, we can predict the chance that a train of a 
certain length appears to be one long null. If we consider that 
we probably miss or misidentify events that start or end with 
pulses in stead of nulls, the occurrence of long apparent nulls 
agrees quite well with the prediction. For groups of length 5, 
for example, the first and last positions are nulls by definition, 
but the middle 3 positions should be nulls, the probability of 
which is (-j-(^^)"'=0.09. This is in good agreement with the 
observed ratio of 1 in 9 that follows from table[2 

This means that at the pulsar the intermittent pulse trains 
and the long nulls are the same. When the short nulls in such 
a train all happen to occur when the pulsar faces us, the effect 
only appears different. 

We now know that during the null-trains, the pulsar blinks 
quickly, producing several nulls that are shorter than a pulse 
period. The continuous increase towards shorter nulls in null- 
length histogram fig.|3already hinted at the existence of a large 
number of nulls shorter than 1 pulse period. Some of the starts 
and ends of these nulls must occur within the pulse window. 
Because of the variability normally already present in the pulse 
shapes, this dying or growing of the flux is not directly rec- 
ognizable, but statistically it results in a lower energy for the 
pulse. In the histogram of pulse energies these 'half-nulls', or 
'half-pulses' should be located in between the energy-values of 
the real nulls and the real pulses (the vertical lines in fig.|2}. 

In the histogram, we do indeed see a number of these low- 
energy pulses (about 0.3%, with energies varying from 10 to 
40% of the average pulse energy). For a pulse to have 10-40% 
of the normal energy, the null should start (stop) in the last 
(first) 20% of the pulse window, so within + 2°. If the positions 
of the starts and stops of null are random, we would expect 
to see this happen in ^=0.5% of pulses, in good agreement 
with the values found above. If these low-energy pulses are re- 



ally half-nulls, we would expect them to occur preferentially 
near null trains, as normal nulls do. We therefore calculated 
the mean distance to a null for both normal and low-energy 
pulses. We find that low-energy pulses are on average about 
30% closer to a null than normal pulses, so from this we con- 
clude that at least some of the low-energy pulses are actually 
half-nulls, where we see the pulsar turning on or off while fac- 
ing us. 
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Table 1. Length, pattern and occurrence of all null-groups in 
our sample. By definition, events always start and end with 
nulls. N denotes a null, P a pulse. Symmetric patterns are 
treated as one: NPNN for example is a combination of NPNN 
and NNPN. Within a group, observing a pulse is 1 .2 times more 
likely than observing a null. 
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Profiles around nulls 

The next step in investigating how the pulsar turns on and off 
is to investigate changes in pulse type or shape before and after 
nulls. First, we compare the number of subpulses in normal 
pulses with the number observed in pulses that lie in between 
nulls (cf. pulses 25, 26, 143 and 145 in fig. [T). We find that 
on average the former contain 1.79+0.01 subpulses, the latter 
1 .80+0.05. To compare pulse shapes, we have averaged the last 
and first active pulses around all nulls in our sample, as shown 
in fig. 13 

The profile of the first pulse after the null is roughly as ex- 
pected; as the pulsar sometimes turns back on only halfway 
through the pulse window, we expect the leading edge of the 
profile to diminish compared to the trailing edge, and this we 
indeed see. For the last pulse before the null one would expect 
the opposite (the leading edge should become relatively more 
important as the pulsar sometimes already turns ofi" near the 
end of the null), but this we do not see. The last pulse before 
the null actually shows a stronger trailing part, contrary to what 
we would expect. 

Finally we note that the first and last active pulse more 
clearly show two-peaked profiles than the average profile does. 
As the separation between the peaks is roughly equal to the 
subpulse separation P2, they could be explained by the pulsar 
turning on and off at certain preferred subpulse-carousel po- 
sitions. In this case the subpulses would not smear out in the 
average, as in the normal profile, but add up. We have not been 
able to verify this suggestion in an independent way however, 
so we put it forward with some reservation. 



Profiles around low-energy pulses 

As some of the low-energy pulses are probably partial nulls, we 
also investigated the profiles of these pulses and the pulses sur- 
rounding them. Similarities between the profiles around nulls 
and low-energy pulses, would again suggest they are related. 



We averaged the profiles in and around 80 low-energy pulses, 
as shown in in fig.|6l 

First thing to note is the two-peaked profile of the average 
low-energy pulses themselves, with the peak separation equal 
to the subpulse separation. This indicates that certain subpulse- 
positions are oveiTepresented, a conclusion strengthened by the 
fact that the last and first active pulses are offset by the same 
amount as would be expected based on normal subpulse drift. 
Apparently, normal pulses with certain subpulse positions are 
less bright than average and contaminate the equally dim half- 
null sample. 

Position jumps over nulls 

Looking at the subpulse drift in PSR B08 18-13 in fig.[T] the 
drift direction seems obvious (slowly towards earlier aiTival) 
but in reality it might be different. The subpulses lack clear, rec- 
ognizable traits (width, height or spacing) so instead of mov- 
ing slowly towards earlier arrival, the subpulses might actually 
move more quickly towards later arrival, and still produce the 
same image as fig.fD As we only observe the subpulse position 
once per pulse period, we undersample the subpulse movement 
and cannot determine the unde rlying speeds immediately (see 
also fig. 2 in lvan Leeuwen et alJ2003l) . 

By focusing on drift behaviour around nulls, we hope to 
obtain more information about the subpulse-drift alias mode 
and the rotation speed of the un derlying carousel, like we pre- 
viously did for PSR B0809+74 jvan Leeuwen et alJ2003l) . For 
each null in our sample, we determine the positions of the pre- 
ceding and following subpulses and calculate the phase jump 
of the subpulse system over each null. 

In fig.Et we plot the longitude jump of the subpulse system 
versus null length, for all nulls in our sample. As the subpulses 
are indistinguishable, it is not clear which pre-null subpulse 
corresponds to which post-null subpulse. Therefore, in fig.Et^ 
all the data is mapped into one region that is P2 (the average 
subpulse separation) wide. To visualize the data better, we have 
recreated the subpulse pattern by duplicating the data-points 
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Fig. 5. Average profile around nulls. Black line: average profile 
for normal pulses. Dotted line: average profile for all pulses 
directly before a null. Dashed line: average profile for all pulses 
directly after a null. We have used all (500+) nulls, including 
those in groups. 
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Fig. 6. Average profile around lows, as in fig. [S] Again, the 
black line represents the average pulse profile. The dotted and 
dashed lines are the profiles before and after the low in this 
case. The dashed-dotted line shows the average profile for the 
low-energy pulses. 
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Fig. 7. Subpulse-pattern longitude jumps over continuous nulls. Each mark is the jump over one null. 

a) Mapped into a region from -|P2 to |P2 around the expected location in the continuous-drift case, b) The same data, duplicated 
to show subpulse grouping, with simulated pre-null subpulse positions added for clarity. The vertical dashed-dotted line predicts 
subpulse positions for a total stop of all drift during the null. The dashed lines are the prediction if drifting continuous during the 
null as it did before, for non-aliased drift. The dotted lines show which subpulses one spark makes if the drift is aliased, c) The 
same data, now with the predicted aliased subpulse path (dotted line), and the actual lower driftrate (solid line). 



with offsets of NXP2 in fig.Eb- The real grouping of the sub- 
pulses is now much more apparent. This diagram is analogous 
to fig. [0 but is now thought to show the drifting of the sub- 
pulses during the null, when they are not actually visible. We 
have also added some simulated pre-null subpulse positions in 
the lower half of the diagram, for clarity. 

If, during the null, the subpulse would stop drifting, the sub- 
pulses should follow the vertical dashed-dotted line in fig.Et". 
If the subpulses continue to drift during the null, they should 
follow the dashed line. Neither of the two is a particularly good 
fit. After the null, the subpulses appear to outrun the prediction- 
line, which would indicate the drift-speed has increased during 
the null. 

The subpulse drift alias mode 

For the dashed line in fig.Eb we have assumed that the drifting 
is not aliased, i.e. that one driftband is made by only one spark. 
Another option woul d be that the carousel is ro tating the other 
way (as in fig. 2b in Ivan Leeuwen et alJl2003l) . In this case, a 
single spark does not move towards earlier arrival and follow 
the dashed line in fig. 05. It moves towards later arrival and 
follows one of the dotted lines. 

Compared to these dotted lines of expected positions, the 
post-null subpulses are lagging, which means the drifting slows 
down during a null. This would be consistent with the be- 
haviour of the only other well-studied case of drifting-nulling 
interaction, in PSR B0809-I-74. There the spark-carousel un- 



ambiguously decelerates during a null, even to the point of 
a total stop. If we assume that our two-pulsar sample be- 
haves similarly instead of oppositely, the spark carousels in 
both should slow down during nulls, and therefore the real 
drift direction in PSR BOS 18-13 must be towards later ar- 
rival. It could do so as illustrated by the dotted lines in fig. 
05, but faster too. If we assum e the carousel consists of 20 
sparks, similar to PS R B0943-Hl0lDeshDande & RankinI lll999l) 
and PSR B0809-H74 Jvan Leeuwen et alJl2003l) . we can place 
an upper limit on the carousel rotation time. In the slowest 
aliased case, a subpulse takes 1.3 pulse periods P to reach 
the position of its neighbour (i.e. the vertical spacing of the 
dotted lines in fig. Eb). The whole carousel then rotates in 
1.3 xPxNspark, ~30s. TWs already makes PSR B0818-13 the 
pulsar with the fastest carousel rotation time known; if the sub- 
pulse drift is more strongly aliased, the rotation time is even 
smaller. 

Subpulse drift during intermittent nulls 

If the carousel in PSR B08 18-13 slows down during a null, 
then why doesn't it stop like the carousel in PSR B0809H-74 
does? Again, the answer is in the intermittent nature of the ap- 
parent long nulls. What appears to be a long null is actually 
a series of short nulls: during each of these the carousel may 
slow down, only to speed up again in between. On average, 
this comes down to a steady lower carousel speed during nulls. 
In fig.El this would translate to a straight line through the post- 
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null subpulse groups, and this is in good agreement with the 
data, as the full line in fig. 0; shows. If we assume the alias 
mode as illustrated by the dotted line in fig.^t, the average drif- 
trate (full line) during a long apparent null is about 0.87+0.05 
of the normal driftrate (dotted line). Using the pulse/null ratio 
of 1.2 we found previously, this means that during the short 
nulls that compose this long apparent null, the drift must be 
less than half the normal driftspeed to produce the apparent 
driftspeed observed. 

If this is correct, and the slow drift-rate seen in the apparent 
long nulls is due to the subpulse drift starting and stopping ev- 
ery time a short null changes into a short burst, we must find the 
same lower driftrate over the intermittent null-groups. In fig.lH 
we show the subpulse pattern over these intermittent nulls in 
the same fashion as fig. |S] did for continuous nulls. Again we 
note that the subpulses arrive earlier than expected. The simi- 
larity in figures0and|S]is striking; even in their drift-rate, the 
continuous and intermittent nulls behave identically. 

4. Discussion 

In a rotating-carousel model, the groups of subsequent short 
nulls and pulses can be explained in two ways. In the one we 
have thus far focussed on, the entire carousel quickly stops and 
starts emitting. In this model all subbeams pause at the same 
time. In a different model, only some subbeams in a carousel 
stop emitting while the carousel rotation continues, creating 
pulse-null patterns similar to those observed. When the off-part 
of the carousel rotates into view, one observes a null. One pulse 
period later, an on-part might have moved into the pulse win- 
dow to produce a normal pulse. Such a configuration would 
have to meet some very specific criteria, though. If the pattern 
is caused by on and off-parts of the carousel turning in and out 
of view, one should also observe half-on half-off" pulses, with 
fewer subbeams. Yet we have found that the pulses in pulse- 



null groups contain exactly as many subpulses as normal pulses 
do. Only a carousel with many subbeams (roughly several hun- 
dreds) at a very high rotation speed could begin to explain this 
lack of observed on- to off'-part transitions, but with a subpulse 
separation of 4°, the inferred size of the total carousel appears 
unrealistic. 

When we compare pulsars B0809-H74 and B0818-13, the 
similarities are striking. Both shine normally for 98-99% of 
the time, to be occasionally interrupted for several seconds and 
then resume their normal emission. Both show similar changes 
in driftrate during nulls. It seems clear we see the same mech- 
anism at work in these two pulsars. Yet for PSR B0809+74 
we know the long nulls must be real, as t he subpulse drift 
stops during these nulls Jvan Leeuwen et al. 200 3). while in 
PSR B08 18-13 we never see long nulls. There we observe pe- 
riods of intermittent emission and quietude, in groups about as 
long as the nulls in PSR B0809+74. 

Apparently, the mechanism that causes the emission in 
PSR B0809+74 to fully turn off, only reaches a border-line 
state in PSR B0818-13. Why the chances for null or pulse oc- 
currence within such a state should be almost equal is unclear. 
The long timescales involved (even in the intermittent nulls) 
remains the biggest problem. In a polar-gap model (the basic 
physics of which are already debatable) for example, the po- 
tential gap over the magnetic pole can build up or discharge in 
microseconds (Ruderman & Sutherland 1975) . Why the pulsar 
would blink on a timescale of seconds, for 10 seconds in total, 
as it does in PSR B08 1 8 - 1 3 remains unclear; we have not been 
able to explain this behaviour with any current pulsar theory. 
What one can do, is study the behaviour of different pulsars to 
find underlying patterns, which is what we have now done for 
PSR B0809-H74 and PSR B0818-13. 

The carousel rotation times found in PSR B0943-I-10 (41 
seconds) and PSR B0809-I-74 (> 200 seconds) pose a similar 
time-scale problem; in a simple polar-gap theory, the rotation 
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times are roughly several seconds (iRuderman & Sutherland 
1 19751) and the rotation times found thus-far are much larger. 
The upper limit of 30 seconds we find for PSR BOS 18- 13 ex- 
tends the sample, back towards the previously expected values. 
A theory like the abovementioned also predicts that the rotation 
times scale with pulse period P and magnetic field strength B 
as ^, which implies that the rotation time of PSR BOS 18- 13 
should be 3 times longer than that of PSR B0S09+74, instead 
of the more than 10 times shorter we find. Our determination 
of the carousel rotation time in PSR BOS IS- 13 is therefore a 
step to finding the dependencies of carousel rotation times on 
other pulsar parameters empirically. 

5. Conclusions 

Both the apparently long nulls and the groups of intermit- 
tent nulls can be quantitatively explained by trains of short 
nulls and pulses. The subpulse-drift behaviour over the two 
is identical: during the nulls, the subpulse drift appears to 
speed up. Assuming the underlying sparks slow down, as they 
do in PSR B0S09-I-74, we conclude that the subpulse drift 
in PSR BOS IS- 13 is aliased. This then leads to a carousel- 
rotation time of less than 30 seconds, the fastest one found to 
date. 
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